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Abstract 



The classical theory of Laplace is not suitable for describing the behavior of micro- 
scopic bubbles. The theory of second gradient fluids (which are able to exert shear 
stresses in equilibrium conditions) allows us to obtain a new expression for surface 
tension and radius of these bubbles in terms of functionals of the chemical poten- 
^r) . tial. This relationship allows us to generalize the results of Cahn-Hilliard (1959) and 

Tolman (1948). 
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1 Introduction 



The microscopic bubbles we consider in this paper contain enough molecules 
to be modelled as a continuous strongly non-homogeneous system. To the cias- 
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sical expression for the free energy valid for homogeneous continua has to be 
added a term depending on the variation of mass density. Usually the addi- 
tional term is assumed to depend on the gradient of density only (Cahn, 1959; 
Serrin, 1986). Hence the expression for stress tensor differs from that valid for 
elastic fluids as even in equilibrium conditions shear components also appear. 
The form of the balance equation of force valid for second gradient fluids is 
found in Germain (1973), Casal (1962). In 1948, Tolman established, by using 
a Gibbs-like approach to Laplace theory, some relationships between the ra- 
dius and the surface tension of bubbles in equilibrium with their liquid phase. 
This work produced many theoritical and experimental investigations concern- 
ing bubbles of quasi- molecular dimensions (Fisher, 1980; Kumar, 1991). In this 
case the theory of Laplace is not suitable as the vapour bubble consists mainly 
or exclusively of the interfacial phase. However, to interpret the experimental 
evidence it is still important to attribute a radius and an energy to micro- 
scopic germs: we do this by means of an equivalent model of Laplace type. 
We note that Cahn and Hilliard did not consider the mechanical aspect of the 
nucleation phenomenon. They only considered the thermodynamic pressure 
(i.e. the spherical part of stress tensor deriving from the classical expression 
for free energy) while in Casal (1985) is shown, for second gradient fluids, the 
existence of a capillary non-spherical stress tensor whose trace includes - but 
does not reduce to - quoted thermodynamic pressure. The analysis of some 
preliminary numerical calculations allows us to conclude that the theory of sec- 
ond gradient fluids, reinterpreted through a comparison with Laplace theory, 
leads to predictions closer to the experimental evidence than those available 
in the literature for radii close to the critical one. A further improvement is 
conceivable by considering non-constant capillarity coefficient as suggested by 
De Gennes (1981). 



2 Equilibrium of bubbles and Gibbs rule by the second gradient 



The balance equations of force for capillary fluids are available in the literature 
(see e.g. Serrin, 1986). The most simple model taking into account the non- 
homogeneity effects introduces a unique additional physical constant C (Casal 
& Gouin, 1985). Its value in mks units is very small: its effects are notable only 
inside interfaces. Balance of forces, in equilibrium conditions, reads (r, p, Q 
denoting stress tensor of non-dissipative capillary fluids, mass density and 
potential of body forces): 



theory 



div t — p grad f2 = 0, 
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with 



-p Id + C grad p ® grad p, where p = V — C h Cp Ap, 



with denoting by V and A the thermodynamic pressure and the Laplace op- 
erator respectively. Equation (1) implies: 

gradP + pgrad(fi-C Ap) = 0. (2) 

With assuming isothermal equilibrium conditions, relation (2) becomes: 

p(p)-C Ap + n = C ste , (3) 

where p is the specific free enthalpy (i.e. the chemical Gibbs potential) relative 
to the fluid in a homogeneous state with mass density p and temperature T. 
If fl is negligible, the equilibrium of a bubble surrounded by its liquid phase 
with density p { is represented by a spherically symmetric profile of density 
satisfying (Rocard, 1952; Cahn & Hilliard,1959): 

Equation(4) has to be supplemented with appropriate boundary conditions: 

- (a) because of spherical symmetry, the derivative of p vanishes at the origin, 

- (b) as we have assumed that the bubble is surrounded by a homogeneous 
liquid, the derivative of p vanishes at infinity. 

If pf is the density of the liquid phase in the equilibrium state with plane inter- 
face and the function p is smooth enough, then the theory of Fuchs equations 
(Valiron) implies that: 

(i) for every p { in an interval p m , pf , Eq. (4) and conditions (a) -(h) uniquely 
determine an increasing mass density profile p(r) and in particular its value 
p v at the origin: we say that p(r) satisfies the capillary fluids version of the 
Gibbs Phase Rule] 

(ii) the function p is twice differentiable at the origin; 

(iii) as -j- Pi > 0, p converges at least exponentially to p l when r tends to 

dp 

infinity. With denoting with a prime the radial derivative and by tp the free 
energy per volume, Eq. (4) implies: 

9C 

Cpi'p* + ^p' 2 = p{p)p' - p( Pl ) P ', (5) 
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and by integrating 

C r 2C 

izP' 2 + / — P'\x) dx = iP(p) - iP(p v ) + pi( Pl )(p v - p). (6) 
2 Jo x 



3 Nucleation energy of bubbles 



We deal here with bubbles which are small with respect to a typical size of 
the liquid phase. More precisely, a mass density field p, whose mean value 
is p Q , in a domain V of volume v(V) represents a small bubble if Pl and 

P{%) ~ Pi 



e C 1 exist such that B = {x G V / 



note that, as / (p — p ) dV = 0, 
Jv 



Po 



> e\ satisfies < £■ We 

v(V) 



Po- Pi 



Po 



= 0(e). We now evaluate (in the 



absence of capillarity effects and at the first order of approximation in e) the 
difference between the free energy of such a small bubble configuration and 
the homogeneous configuration of density p . This difference, which is part of 

the nucleation energy of the bubble, is w = I (ip(p) — fp(Po)) dV. 

Jv 

As / p dV = / p dV, the energy w = / [ip(p) — ap— ip(p ) + ap ] dV does not 

Jv Jv Jv 

depend on the choice of the constant a. The total energy is well defined but 

its localization is somewhat arbitrary. This indetermination is related to the 

arbitrariness in the choice of chemical potentials. With the choice a = p(pi) 

the energy is localized inside B. Indeed, with using the previous estimation 

for p — p l we obtain: 

w = [ m P ) - p( Pl )(p - Pl ) - i>{ Pl )\ dV + 0(e 2 ). 
Jb 

To the non-capillary part of nucleation energy of a bubble in an unbounded 
domain we have to add the interfacial (capillary) energy. In the Gibbs-Laplace 
theory we have for bubbles of radius R and surface tension a: 



W = 47TR 2 <T + -7Ti? 3 [^(pJ - + P(Pl)(Pl ~ Pv)\- 



(7) 



In the theory of second gradient fluids, following Cahn and Hilliard, we have: 

W=f m P ) - HPi) - P(Pi)(P - Pi) + ^(gradp) 2 ] dV. (8) 
Jv I 

Let V =pp — ifi denote the thermodynamic pressure. The conditions 

p(Pi) = P(pv), (9) 
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npi)-np v )=-% 



(10) 



valid - for isothermal equilibrium - only in Laplace theory, transform Eq.(7) 
into: 

W = -irR 2 a, (11) 
3 

and imply that the nucleation energy of the bubble is the third of the creation 
energy of its interface. We can extend this result to the theory of second 
gradient fluids. Let (f>(p) = ip(p) — ip(pi) — p(pj)(p — pj; by multiplying Eq. 
(6) by r 2 and integrating it over [0,oo] we get: 



r r2 \2 p ' 2 ~ ^ (p) ) dr + r r2 ^ + Jo ~ p/2 ^ dx 

by integrating by parts and using Eq. (6) again, we obtain: 



dr = 0, (12) 



= o. 



(13) 



Because of (i) - (iii) of Sect. 2, the last term vanishes so that we obtain 

roo / C) \ A roo 

W = 4n r 2 -p' 2 + Mp) ) dr=-7r Cr 2 p' 2 dr. 
Jo V 2 / 3 Jo 



(14) 



With denoting by R 2 m the mean value of r 2 with respect to the measure p' 2 dr, 
Eq.(14) reads 

w = t n ri r cpn dr - 

For a large enough bubble (i.e. when p { tends to pf): 
(i) R m is the radius, 



(15) 



(ii) Equation (15) reduces to Equation (11), as surface tension for plane in- 



terface is 



roo 

/ c P ' 2 

Jo 



dr. 



4 Comparison between Laplace and second gradient theories. Equiv- 
alent bubbles 



In second gradient theory the stress tensor in the center of a spherical bub- 
ble takes the value r = —p v Id where p v = V(p v ). As Equation (4) implies 
CpAPv = Pv (p(Pv) ~ P(Pi)), w e have 

Pv-m = 4>(pi) - i>(p v ) + Kpi)(Pv - Pi)- ( 16 ) 

Let us note that this difference is not equal to the corresponding difference 
of thermodynamic pressures as, for microscopic bubbles, p(pi) differs from 
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n{p v ). As experimental results (Fisher, 1980) deal with measures of stresses, 

then we have to use p v —p\ instead of V(p v ) —V(pi) in the comparison between 

Laplace and second gradient theories. We can now define the surface tension 

and the radius of a bubble by identifying the nucleation energies and the 

2a 

pressure differences computed in both theories. Indeed p v — pi = — and 

4 roo 4 
-7r / Cr 2 p' 2 dr = -nR 2 a imply 
3 Jo ' 3 



R 



R 



2C 



r 2 p' 2 dr 



bP(pi) - i>(Pv) + M(pv - pi)\ * 



(17) 



and 



C 
1 



r 2 p' 2 dr 



[Hpi)-^(Pv)+P(Pi)(Pv-Pi)? 



(18) 



5 Conclusion 



Let us notice that our treatment is based on Eq. (1) obtained by a continuum 
model. However, van Kampen (1964) obtained the same differential equation 
by using ideas and methods of statistical mechanics. Eq. (14) has allowed us to 
evaluate the nucleation free energy of a bubble, once the mass density profile is 
known, i.e. once the liquid phase density p l has been assigned. Equation (16) 
allows us to understand the difference between thermodynamic and stress 
pressure. These two equations do not depend on a particular constitutive law 
and determine the radius and surface tension for a microscopic bubble. This 
is done by identifying the pressure jump and nucleation energy of Laplace 
theory with the corresponding quantities in the second gradient theory. This 
method has never been used in literature: for example, Cahn (1959) and Evans 
(1979) used for surface tension the expression valid for plane interfaces, and did 
not investigate the relationship between surface tension and radius. Our Eqs. 
(17)-(18) give such a relationship in an implicit form: in fact, both surface 
tension and radius depend on the density profile, which, in turn, depends 
on p t . When the radius tends to infinity, our expression for surface tension 
reduces to that for plane interfaces. We observe that the relationship between 
a and R depends on thermodynamic potential. Numerical calculations were 
performed by Dell'Isola, Gouin and Rotoli (1996) using equations (17)-(18) 
and van der Waals-type potential proposed by van Kampen (1964), Rocard 
(1967) and Peng (1976). It was able to observe that, for bubbles whose radius 
is close to the critical one, the predicted variations of the surface tension (for 
water and cyclohexane at 20°C) are in good agreement with the experimental 
observations (Katz et al 1976; Fisher & Israelachvili, 1980). This is not the 
case for the theoretical treatment proposed by Tolman (1947) and Kumar et 
al (1991). 
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Abridged French version 

Rayon et tension superficielle des bulles microscopiques 
en theorie du second gradient 

Resume 

La theorie de Laplace est inadaptee pour etudier les bulles de dimensions moleculaires 
La theorie des fluides doues de capillarite interne nous permet de proposer une ex- 
pression de la tension superficielle et du rayon des bulles comme fonctionelles du po- 
tentiel chimique. Cette expression, en accord avec l'experience, ameliore les resultas 
obtenus par Cahn-Hilliard et Tolman. 

Version francaise abregee 

Pour etudier le comportement d'un fluide fortement heterogene on ajoute a l'expres- 
sion de la densite d'energie libre d'un fluide homogene un terme limite, le plus sou- 
vent, a un developpement au second ordre en gradients (Cahn, 1959; Serrin, 1986). 
Les equations du mouvement de ces fluides sont obtenues par la theorie du second 
gradient (Casal, 1961, Germain, 1973). 

Tolman (1948) a etudie la tension superficielle des bulles en fonction de leur rayon. 
Des mesures experiment ales ont ete faites pour des bulles de dimensions quasi mole- 
culaires (Fisher, 1950; Kumar 1991). Pour de telles dimensions la theorie de Laplace 
est disqualified car la bulle de vapeur devient l'interface elle-meme. II est neanmoins 
important de donner une dimension et une energie pour les germes microscopiques 
puis de les comparer a un modele equivalent de type Laplace. On peut ainsi in- 
terpreter les resultats des experiences. Le modele le plus simple de fluide du sec- 
ond gradient ne fait intervenir qu'une seule constante physique supplementaire C. 
L'equation de l'equilibre est alors l'equation (1) (Serrin, 1986) ou r s' interprete 
comme le tenseur des containtes, Q,p et P designent respectivement le potentiel 
des forces de masse, la densite de masse et la pression hydrodynamique. Si les forces 
de masse sont negligeables, l'equilibre isotherme d'une bulle supposee spherique 
dans une phase liquide de densite p t est represente par une solution de l'equation 
(4). Pour chaque valeur p l dans une plage convenable on obtient un profil croissant 
unique et une valeur de la densite p v au centre de la bulle. Le regie de Gibbs est 
done bien verifiee. Nous comparons l'energie libre d'un fluide homogene de densite 
p Q contenu dans un domaine D et l'energie libre de la meme masse de fluide dans 
D formee par un liquide homogene contenant une petite bulle V . Cela nous conduit 
aux expressions (7) et (8) pour l'energie de nucleation W d'une bulle dans un do- 
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maine infini, respectivement en theorie de Laplace et en theorie du second gradient. 
Les conditions d'equilibre liquide-vapeur en theorie de Laplace permettent d'ecrire 
W sous la forme de l'equation (11). La serie d'egalites (5), (6), (12), (13), (14) mon- 
tre que ce resultat se generalise sous la forme de l'equation (15) ou designe la 
valeur moyenne de r 2 relativement a la mesure p' 2 . En theorie du second gradient, 
le tenseur des contraintes est spherique a l'origine ainsi qu'a l'infini mais, alors qu'a 
l'infini la pression se confond avec la pression thermodynamique pi = P(pi), cela est 
faux a l'origine ou p v = P(p v ) — Cp v Ap v . En identifiant la difference de pression 
entre l'origine et l'infini ainsi que l'energie de nucleation dans les deux theories, on 
obtient les definitions (17) et (18) pour le rayon R et la tension superficielle a de 
la bulle. Notre raisonnement est base sur l'equation (1) obtenue par un modele de 
mecanique des milieux continus. II faut noter que van Kampen (1964) a obtenu la 
meme equation par des considerations de mecanique statistique. L'equation (14) 
nous permet d'evaluer l'energie de nucleation a partir du profil de densite done de 
Pi- La determination du rayon et de la tension de surface d'une bulle microscopique 
par identification de l'energie de nucleation et du saut de pression donne par la 
theorie de Laplace avec les quantites correspondantes en theorie du second gradient 
n'a jamais ete utilisee dans la litterature: par exemple Cahn (1959) et Evans (1979) 
ont utilise une expression de la tension superficielle valable pour une interface plane 
uniquement et n'ont pas etudie la relation liant le rayon et la tension superficielle. 
Les equations (17) et (18) permettent d'obtenir une telle relation sous forme im- 
plicite: le rayon et la tension superficielle dependent tous deux du profil de densite 
qui lui-meme depend de p t . Notre expression de la tension de surface se reduit, quand 
le rayon de la bulle R tend vers l'infini (i.e. quand p l tend vers la valeur pf valable 
pour une interface plane), a l'expression usuelle (Cahn, 1959). L 'etude des variations 
de la tension superficielle en fonction du rayon necessite le choix d'une loi d'etat. 
Dell'Isc-la, Gouin et Rotoli (1996) ont utilise les potentiels chimiques du type de 
van der Waals proposes par van Kampen (1964), Rocard (1967) ou Peng (1976), 
pour effectuer des calculs numeriques: ils montrent que, pour des bulles de rayon 
voisin du rayon critique, la variation prevue de la tension superficielle correspond 
aux mesures experiment ales de Katz (1976) et Fisher (1980) alors que les previsions 
de Tolman (1947) et Kumar (1991) ne sont pas realisees. 
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